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The recent discovery of topological insulators with ex- 
otic metallic surface states has garnered great interest 
in the fields of condensed matter physics and materi- 
als science. 1 A number of spectacular quantum phenom- 
ena have been predicted when the surface states are un- 
der the influence of magnetism and superconductivity}^ 
which could open up new opportunities for technological 
applications in spintronics and quantum computing. To 
achieve this goal, material realization of topological insu- 
lators with desired physical properties is of crucial im- 
portance. Based on first-principles calculations, here we 
show that a large number of ternary chalcopyrite com- 
pounds of composition I-III-VI 2 and II-IV-V 2 can real- 
ize the topological insulating phase in their native states. 
The crystal structure of chalcopyrites is derived from the 
frequently used zinc-blende structure, and many of them 
possess a close lattice match to important mainstream 
semiconductors, which is essential for a smooth integra- 
tion into current semiconductor technology. The diverse 
optical, electrical and structural properties of chalcopy- 
rite semiconductors, 6 and particularly their ability to host 
room-temperature ferromagnetism, 7 8 make them appeal- 
ing candidates for novel spintronic devices. 

Topological insulators are new states of quantum matter 
characterized by the so-called Z 2 topological invariant as- 
sociated with the bulk electronic structure. 9 11 Their band 
structures is usually characterized by a certain type of band- 
inversion that involves the switching of bands of opposite par- 
ity around the Fermi levelJSS Early material realizations of 
to polog ical insulators are limited, including the Bii^Sbz al- 
loyEHniand the binary tetradymite semiconductorsP^KlMoti- 
vated by their application potentials, the search for new topo- 
logical insulators has moved into ternary compounds as they 
provide greater material flexibility. 18 24 For example, it was 
proposed that in the Heusler family 18 20 the topological insu- 
lator can be made with coexisting magnetism, a much desired 
property for spintronic applications. However, a clear short- 
coming of the Heusler family is that in their native states these 
materials are either zero-gap semiconductors or semimetals. 
Although alloying or proper strain can be used to achieve the 
insulating behavior, this creates extra complexity in material 
growth and could introduce detrimental effects upon doping. 
In contrast, the chalcopyrite compounds reported in this work 
have an intrinsic band gap due to its crystal structure, there- 
fore is free from these problems. 

The ternary ABCn chalcopyrite compounds of composition 
I-III-VI2 or II-IV-V2 are isoelectronic analogs of the II- VI or 
III-V binary semiconductors, respectively. 6 The crystal lattice 
of chalcopyrites is described by the space group (/42<f), 
which can be regarded as a superlattice of the zinc -blende 
structure with small structural distortions. As illustrated in 



Fig. 1, the chalcopyrite unit cell is essentially formed by two 
cubic zinc-blende unit cells, in which the A and B atoms are 
ordered on the two different cation sites. In view of the over- 
all structural similarity between the ternary chaclopyrites and 
their binary analogs, the electronic states of the former are 
expected to closely resemble those of the latter. In particu- 
lar, we anticipate that some of the chalcopyrites might fall 
into the same topological class as the topologically nontriv- 
ial 3D-HgTep^^l which is a II- VI zinc-blende semiconductor. 
Furthermore, the cubic symmetry in chalcopyrite compounds 
is explicitly broken due to the AB cation ordering, accompa- 
nied by two additional structural modifications (Fig. lb): the 
tetragonal distortion characterized by the ratio of the lattice 
constants rj = c/2a, and the internal displacement Su of the 
anions away from their ideal zinc -blende sites. In an ideal 
structure, r\ = 1 and Su = 0. As we show below, this inherent 
symmetry reduction plays the same role as the uniaxial strain 
required in cubic topological materials such as 3D-Hg TePand 
the Heusler familyP^HImaking it possible to realize the topo- 
logically insulating phase in the native state of chalcopyrite 
compounds. 

To illustrate the band topology of the ternary chalcopyrites, 
we begin with their binary analog, HgTe, and follow its band 
structure under an adiabatic transformation during which the 
zinc -blende structure is gradually evolved into the chalcopy- 
rite structure (Fig.|2^-c). The first-principles relativistic band 
structure is obtained by using the full-potential linearized 
augmented plane-wave method implemented in the WIEN2K 
package. 25 The exchange-correlation potential is treated using 
the modified Becke-Johnson potential, 26 which is a semi-local 




FIG. 1: Comparison of the zinc-blende and the chalcopyrite 
structure, a, The zinc-blende compound HgTe. b, The chalcopyrite 
compound ABC-z. The internal displacement of the anion is defined 
as 5u = (Rac ~ R%c)/ a2 > w h ere Rac and Rbc are the bond 
lengths between the anion C and its two nearest A and B cations. 
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FIG. 2: Band structures of HgTe in comparison with AuTlTe2 and AuTlS2. a-c, Band structures of HgTe (a), HgTe in an ideal chalcopyrite 
supercell (b), and HgTe in a chalcopyrite cell with 77 = 1.016 and Su = —0.018 (c). The size of red dots denotes the probability of s-orbital 
occupation on the cations. The insert in c shows the energy gap opened by chalcopyrite distortion. d,e, Band structures of AuTlTe2 (d) and 
AUTIS2 (e). Both display an inverted band order similar to HgTe, indicating a nontrivial band topology. AE in e indicates the band inversion 
strength. 



potential designed to accurately predict the band gap.^The 
band structure of HgTe features a distinctive inverted band or- 
der at the r point, where the s-like Tq states sit below the 
four-fold degenerate p-like T 8 states (Fig. [2^). Away from the 
r point, the valence and conduction bands are well separated 
without crossing each other. Since the band inversion occurs 
only once throughout the Brillouin zone, HgTe is in a topolog- 
ically nontrivial phase. 13 The four- fold degeneracy of the Tg 
states is protected by the cubic symmetry of the zinc -blende 
structure, thus in its native state HgTe is a zero-gap semicon- 
ductor. The band structure calculation is then repeated for 
HgTe in an ideal chalcopyrite supercell (Fig. |2]3), which is 
four times larger than the zinc -blende one. During this pro- 
cedure, every set of four different wavevectors of the original 
zinc -blende Brillouin zone fold into a single point of the chal- 
copyrite Brillouin zone. In particular, states at the original V 
point fold back to themselves, therefore the low-energy elec- 
tronic properties are still dominated by states around the T 
point. In the final step, to mimic the cubic-symmetry break- 
ing in chalcopyrites, we let the system undergo a tetragonal 
distortion plus an internal displacement of the two Te atoms. 
This structural change lifts the degeneracy of the Tg states and 
opens up an energy gap, while leaving the inverted band order 
intact (Fig. |2};). The situation here thus is exactly like HgTe 
under a uniaxial strain,LU13U an( j me resu hing compound, an 
artificial chalcopyrite Hg2Te2, is a topological insulator. 

The above analysis of the tetragonally distorted HgTe sug- 
gests that the same Z2 band topology can be realized in real 
ternary chalcopyrites. This is confirmed by our extensive first- 
principles band structure calculations over this family of ma- 
terials. There are two possible scenarios, with the correspond- 



ing band structures of representative systems shown in Fig.|2jl 
and 2e. The most natural generalization of HgTe to chalcopy- 
rites is the ternary compound AuTlTe2, which can be obtained 
by replacing the two Hg atoms in Hg 2 Te 2 with an isoelec- 
tronic pair of Au and Tl atoms. Conceptually, this process can 
be described by an adiabatic transformation in which the nu- 
clear charges on the two Hg sites are continuously changed to 
that of Au and Tl, respectively. However, such a substitution 
causes the conduction bands to drop below the Fermi level at 
the N point, making the system effectively a semimetal. Note 
that a Z 2 topological invariant can still be defined for the va- 
lence bands in AuTlTe2 because they are separated from the 
conduction bands by local energy gap throughout the Bril- 
louin zone. Our further study shows that the semimetallic 
behavior can be removed by choosing a lighter element S 
in place of Te — we find that a global band gap of 0.14 eV 
can be realized in AUTIS2 (Fig. |2£). Remarkably, this band 
gap is much larger than typical band gaps opened by uniax- 
ial strain in Heusler compoundsJi^2l an( j j s a j so we jj arJ ove 
the energy scale of room temperature. In both AuTlTe2 and 
AUTIS2, the inverted band order at the T point is retained, 
with no band inversion at other places. We therefore conclude 
that, in their native states, AUTIS2 is a topological insulator 
with a sizable band gap, and AuTlTe2 a topological metal. On 
the other hand, compounds with normal band order are always 
semiconductors — our calculation did not find any semimetal- 
lic chalcopyrite with normal band order. 

So far our argument has been based on the isoelectronic 
relation between the ternary chalcopyrites and their binary 
analogs; it is much desired to give a definite proof that the 
band topology of AuTlTe2 and AUTIS2 is Z2 nontrivial. Be- 
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FIG. 3: Inverted band strength for various chalcopyrites as a function of the lattice constant. Open symbols mark the compounds whose 
lattice constant was reported in the literature. For the rest of the compounds their equilibrium lattice constants are obtained by first-principles 
total energy minimization. When AE < 0, squares mark the compounds that are topological insulators and diamonds mark the systems that 
are topological metals. Shaded areas indicate materials that are expected to be closely (±2%) lattice matched to either GaAs, InAs, or InSb. 



cause of the lack of inversion symmetry in these materials, 
the parity criterion developed by Fu and Kane 13 does not ap- 
ply. Here we directly evaluate the Z 2 invariant in terms of 
the Berry gauge potential and Berry curvature associated with 
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FIG. 4: Band structures of CdSnAs2 without and with a hydro- 
static strain. The application of a hydrostatic strain causes the s-like 
bands (marked by red dots) to jump below the valence band top, 
providing the necessary band inversion that leads to a nontrivial topo- 
logical order. 



the Bloch functions.^ In this approach, the Z 2 invariant is 
interpreted as an obstruction to smoothly defining the wave 
function throughout the Brillouin zone under a gauge that re- 
spects the time-reversal symmetry 121S2J Unlike the parity crite- 
rion, this approach does not require any specific point-group 
symmetry. Following ref. 30, we have carried out the lattice 
computation of the Z 2 invariant based on our first-principles 
method (see ref. 18 and Supplementary Information for de- 
tails). We find that, for all the compounds presented in this 
work, the calculated Z2 invariant agrees with the intuitive 
band-inversion picture. In practice, this translates into an em- 
pirical rule that if the s-orbital originated states are com- 
pletely occupied and below the valence band maximum, the 
compound possesses a nontrivial topological band structure. 

Armed with the above insight, next we explore the band 
topology of the chalcopyrite family of the I-III-VI2 com- 
pounds (I = Cu, Ag, Au; III = In, H; V = S, Se, Te), as well 
as the II-IV-V2 compounds (II = Zn, Cd, Hg; IV = Ge, Sn; 
V = As, Sb). Since the band inversion only occurs at the T 
point, the band order at that point can be used to characterize 
the band topology. In compounds with cubic symmetry, one 
can define the band inversion strength as the energy differ- 
ence between the Tq and Tg states. 19 20 However, in passing 
from the zinc -blende structure to the chalcopyrite structure, 
the original T 8 states split into states with T 7 and T 6 symme- 
try, which typically form the top set of the valence bands and 
the bottom set of the conduction bands. We therefore define 
the band inversion strength AE as the energy difference be- 
tween the s-orbital originated T 6 states and the valence band 
maximum at the T point. Figure|3]shows AE for a number of 
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chalcopyrite compounds as a function of the lattice constant. 
Materials with AE > are topologically trivial and those 
with AE < are either topological insulators or topological 
metals. As shown in Fig. [3] there are quite a few chalcopyrite 
topological insulators with a close lattice matching to several 
mainstream semiconductors. Finally, we note that the impres- 
sive collection of chalcopyrite topological insulators can be 
further extended by proper strain engineering. Similar to the 
Heusler family ] 18 " 20 ' the topological phase is sensitive to the 
lattice constant, which offers certain degree of tunability of 
these materials. For example, a 7% change in the lattice con- 
stant converts the trivial semiconductor CdSnAs2 into a topo- 
logical insulator (Fig. |4j». 

The interplay between topological surface states and ferro- 
magnetism has been under extensive investigation because of 
the possibility of electric control of magnetization with low 
power consumption.^ The wide availability of chalcopyrite 
semiconductors, the control of topological order via lattice 
tuning, and their excellent prospect for room temperature fer- 
romagnetism, make them the ideal platform to further investi- 
gate these novel topological phenomena and could enable the 



next-generation spintronic devices. 
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